Van Hove singularities in the paramagnetic phase of the Hubbard model: a DMFT study 
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Using the dynamical mean-field theory (DMFT) we study the paramagnetic phase of the Hubbard model with 
the density of states (DOS) corresponding to the three-dimensional cubic lattice and the two-dimensional square 
lattice, as well as a DOS with inverse square root singularity. We show that the electron correlations rapidly 
smooth out the square-root van Hove singularities (kinks) in the spectral function for the 3D lattice and that 
the Mott metal-insulator transition (MIT) as well as the magnetic-field-induced MIT differ only little from the 
well-known results for the Bethe lattice. The consequences of the logarithmic singularity in the DOS for the 2D 
lattice are more dramatic. At half filling, the divergence pinned at the Fermi level is not washed out, only its 
integrated weight decreases as the interaction is increased. While the Mott transition is still of the usual kind, 
the magnetic-field-induced MIT falls into a different universality class as there is no field-induced localization 
of quasiparticles. In the case of a power-law singularity in the DOS at the Fermi level, the power-law singularity 
persists in the presence of interaction, albeit with a different exponent, and the effective impurity model in the 
DMFT turns out to be a pseudo-gap Anderson impurity model with a hybridization function which vanishes at 
the Fermi level. The system is then a generalized Fermi liquid. At finite doping, regular Fermi liquid behavior 
is recovered. 

PACS numbers: 71.27.+a, 71.30.+h. 72.15.Qm 



I. INTRODUCTION 

For many materials, it is permissible to consider each elec- 
tron as moving essentially independently in a static periodic 
effective potential which takes into account the interactions 
between the electrons. This point of view received theoreti- 
cal support through the density-functional theory (DFT)i and 
the Kohn-Sham Ansatz-, which consists of replacing the full 
many-body problem with an auxiliary independent-particle 
problem. DPT calculations are remarkably accurate for wide- 
band systems, such as simple metals, many semiconductors 
and insulators, but they are less appropriate for strongly cor- 
related electron systems, such as some transition metals, lan- 
thanides and their compounds^ "^. 

For non-interacting systems, the density of states (DOS) 
counts the number of the single-particle levels that may be 
occupied per unit energy. It can be defined as 



be present in any band structure^. In two-dimensional sys- 
tems, for example, the saddle point in the dispersion gives 
rise to a logarithmic divergence in the DOS. The van Hove 
singularities are thought to be particularly important for the 
physics of low-dimensional systems, for example in oxide 
superconductors^>>2i2i2ii£. 

For strongly-interacting systems, the concept of the single- 
particle levels is not very useful and one should resort to the 
techniques from the many-particle theory. In particular, the 
equivalent of the DOS is the local spectral function 

p{^) = --Im [Gioc(^^ + iS)] , (3) 
n 

where G\oc{i^) is the local (k-averaged) Green's function 

GiocH = ^^GkM, (4) 



(1) 



where k indexes the N single-particle levels with energies cu- 
For an infinite system, the sum goes into an integral 
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Any smooth periodic function must have critical points where 
the gradient vanishes. In a periodic structure such as a crys- 
tal, Ek is periodic in the reciprocal space, therefore po{lu) 
will necessarily have singularities arising from minima, max- 
ima and saddle points of Ck- For topological reasons, a cer- 
tain minimum number of these van Hove singularities must 



with Gk(2) = ((ck; cj^))z the momentum-resolved Green's 
function (electron propagator). In the absence of interactions, 
Im [Gk(i^ + iS)] = — 7r(5(w— ek) andEq. ([3]) reverts to Eq. ([T]). 
The interactions modify the propagation of electrons, thus 
the spectral function p{uj) differs significantly from the non- 
interacting DOS po{oj) and, in particular, any sharp features 
such as van Hove singularities will smooth out due to broad- 
ening (finite life-time) effects. 

The presence of a diverging DOS at the Fermi level en- 
hances instabilities towards various ordered states, in par- 
ticular antiferromagnetism and ferromagnetisniii, as well as 
superconductivity'^. In addition, the underlying paramagnetic 
metal phase itself is expected to have unusual properties. In 
this work, we address this regime using the dynamical mean- 
field theory (DMFT)i^, using the numerical renormalization 
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group (NRG) as the impurity solve r'"^'^^'^^ . We will study the 
spectral function of the Hubbard mode l^^''^''^ in the param- 
agnetic phase on the three-dimensional simple cubic lattice 
(which has square root singularities at the edges of the band 
and two further square root singularities inside the band) and 
on the two-dimensional square lattice (where a logarithmic 
singularity is located in the center of the band). We focus on 
the effects (in particular on metal-insulator transitions'*) where 
local physics plays the key role, thus the use of the DMFT as 
an approximate method to study finite-dimensional systems is 
justified. For completeness, we also consider the case of a 
DOS with an integrable power-law singularity, to wit an in- 
verse square root DOS. While such divergences arise in ID 
problems, we consider here this case using the DMFT purely 
out of academic interest. 



II. METHOD 

In the limit of infinite dimensions or infinite lattice connec- 
tivity, the self-energy in the Hubbard model becomes purely 
local'^'^"'^'. The bulk problem of correlated electrons then 
maps exactly onto a quantum impurity model with a self- 
consistently defined non-interacting bath of conduction elec- 
trons, in this case an Anderson impurity model^^i2i2i2^. To 
solve the effective impurity model, various non-perturbative 
techniques can be used, for example the numerical renormal- 
ization group (NRG) 14, 15.26,27,28.29,30,3 1 which is particularly 
suitable to study the low-temperature limit. The NRG cal- 
culations in this work were performed for the discretization 
parameter A = 2, with the z-averaging over 8 values of the 
twist parameteri^^i^ using a modified discretization scheme 
from Ref. [T^ Spectral functions were computed using the 
density-matrix approach^ and the self-energy trick^^. The 
truncation cutoff was set at i?cutoff = Wlun and the broad- 
ening was performed with parameter a = 0.1. 

The input to a NRG calculation step in the DMFT cycle is 
an effective hybridization function Tcr{u!) which contains full 
information about the coupling between the interacting impu- 
rity site and the effective non-interacting medium. The output, 
as required for the DMFT calculation, is the self-energy func- 
tion I]cr(ijj). The local lattice Green's function is then 



N ^ [iV + Ha - Scr(t^)] - fife 

Po{e-)de 



(5) 

(6) 

(7) 
(8) 



where po{z) is the density of states in the non-interacting 
model, while Go{z) is the associated free-electron propaga- 
tor. The local lattice spectral function is then 
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The self-consistency condition'^ relates the local lattice 
Green's function G\oc,a and the hybridization function as 
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(10) 



The DMFT iteration proceeds until two consecutive solu- 
tions for the local spectral function differ by no more than 
some chosen value. To accelerate the convergence to the self- 
consistency and to stabilize the solutions, one can make use 
of the Broyden mixing^^. 



III. 3D CUBIC LATTICE DOS 

An analytical expression is known for the Green's function 
for the three-dimensional simple cubic latticei^; 
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where K{m) is the complete elliptic integral of the first 
kind. When using these expressions, one should be careful 
to choose the correct convention for the argument of the ellip- 
tic integrals, since several different are in common use. In the 
dimensionless form of Eq. (fTTT i. the band is centered at zero 
and the bandwidth is equal to 6. The behavior of the DOS at 
the band edges is typical for three-dimensional systems: the 
DOS goes to zero as V Alu, where Alu is the distance from the 
band edge. The DOS is thus continuous, however its deriva- 
tive diverges. This feature is shared by the DOS of the Bethe 
lattice with infinite connectivity which is, for this very reason, 
a common model DOS for three dimensional systems. Inside 
the band there are two further square root singularities where 
the DOS is continuous, but the derivatives are discontinuous 
and diverging on one side. 

All four van Hove singularities smooth out as soon as the 
interaction is turned on, see Fig.[T^. At moderate U/W ^ 0.6 
the spectral function already strongly resembles that obtained 
for a structureless DOS such as the one for the Bethe lattice 
with infinite coordination^*^. Nevertheless, the shape of the 
non-interacting DOS does have an effect on the quantitative 
details. The Mott-Hubbard metal-insulator (MIT) transition 
occurs at a value of Uc/W = 1.165, which is to be compared 
with the results for the Bethe lattice, Uc/W = 1.47^^". If the 
result is rescaled in terms of the effective bandwidth defined 
through the second moment of the DOS^" 



Im [Gioc,o-(w -I- iS)] 



(9) 



/ fW/2 

Wcs = 4 J / de e2p(e) « 0.816W, 

V J-W/2 



(12) 
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we obtain Uc/WcS = 1-43, which compares well with the 
Bethe-lattice result. The feature that W^ff is a rather ro- 
bust characteristic preferable to the bandwidth has been noted 
before 
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Figure 1 : (Color online) a) Spectral functions for the Hubbard model 
with the 3D cubic lattice DOS in the paramagnetic phase at half- 
filling, r = 0. Due to the particle-hole symmetry only positive fre- 
quencies are shown. The arrows indicate the direction of increasing 
interaction U. The dashed curve corresponds to a result for the insu- 
lating phase, b) The evolution of the spectral function as a function 
of the magnetic field. The inset shows the field-dependence of the 
magnetization. 



The Hubbard model undergoes a magnetic-field-induced 
MIT in external magnetic field, the main mechanism on the 
Bethe lattice being a field-induced quasiparticle mass en- 
hancement (field-induced localization ) '-^i^^i^^ . The behavior 
for a 3D cubic lattice DOS is very similar: the quasiparti- 
cle peak narrows down and moves slightly away from the 
Fermi level, while the Hubbard bands become increasingly 
spin-polarized and they take the form of the non-interacting 
DOS, however they do not move much, see Fig.[T}5. The MIT 
occurs when the quasiparticle peak vanishes. This happens 
for a field on the scale of the width of the quasiparticle peak 
in the absence of the field. For U /W = 1, for example, we 
obtain g^isBc = 0.099M^, which is to be compared with the 
quasiparticle residue Z = 0.077. 



IV. 2D SQUARE LATTICE DOS 

The free-electron propagator on the two-dimensional 
square lattice is^ 



Go{z) 



-K 



1 



with a z 



Goiz) = 



expansion 

-i(21nz-41n2) 
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(13) 



(14) 



which gives the logarithmic singularity in the density of states 
at the Fermi level for a half-filled system: 



po{uj) 



21ncj - 41n2 



27r2 



(15) 



The divergence at the Fermi level in the DOS is not elimi- 
nated as the interaction is turned on, see Fig. |2l The effect 
of the interaction at low energy scales is to renormalize the 
constant part 41n2/(27r^) to smaller values, while the loga- 
rithmic term keeps the same prefactor We find that the loga- 
rithmic scaling is difficult to achieve numerically to very low 
energies, where spurious features were observed for energies 
below lO^^Pl^. In spite of the diverging spectral function, 
the Mott metal-insulator transition appears to be of the usual 
type. With increasing interaction U, a region of low spectral 
density appears before the onset of the logarithmic peak. We 
find Uc/Wcff — 1.45, in agreement with the standard result (it 
should be noted that for the 2D square lattice DOS we have 
WcS = W). The value of Uc was determined by studying the 
sharp resonances in Iml](aj) in the metal phase which evolve 
into the zero-frequency pole in the self-energy for the insula- 
tor phase^^ . We extracted their position X as a function of the 
interaction U, and solved for Uc in the equation X{Uc) = 0. 

For strong interaction, the Hubbard bands have a strongly 
asymmetric shape with pronounced shoulders near the inner 
band edges, and less pronounced shoulders at the outer band 
edges. These features survive into the insulating phase, which 
is to be contrasted with the behavior in the case of Bethe 
and 3D cubic lattice, where the structure at the inner band 
edge disappears as the transition point is approached^^^. 
The nature of the shoulders observed here is thus different: 
they reflect the discontinuities at the band edges of the non- 
interacting DOS. 

When the non-interacting system is doped, the logarith- 
mic singularity moves away from the Fermi level. Albeit 
the energy resolution of NRG is finite at energies away 
from the Fermi level, it is nevertheless sufficient to study 
the interaction-induced broadening (the energy resolution of 
NRG at finite energies has been recently studied in Ref. 
A plot of the spectral functions as a function of doping is 
shown in Fig. |3] At finite doping, the imaginary part of the 
self-energy is finite at the energy of the impurity level, thus 
the singularity transforms into an asymmetric Lorentzian-like 
peak. The quasiparticle residue (wavefunction renormaliza- 
tion) 
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Figure 2: (Color online) Spectral functions for the Hubbard model 
with the 2D square lattice DOS in the paramagnetic phase at half- 
filling, plotted on a a) linear and b) logarithmic energy scale. The 
arrows show the direction of increasing U . The result plotted with 
dashed curve is already in the paramagnetic insulator phase. 



with increasing field is different from that in the systems with 
finite DOS at the Fermi level; it grows broader rather than 
narrower, while its integrated weight decreases, see Fig. |4] 
As the transition point is approached from below, we find 
that the DMFT calculations no longer converge even when 
the advanced Broyden mixing is used. This signals the com- 
plete absence of solutions to the DMFT equations, rather than 
merely their instability or metastability. This might signal the 
presence of an intermediate phase which must be qualitatively 
different from both the homogeneous partially-spin-polarized 
paramagnetic phase and the fully polarized insulating phase. 




Figure 4: (Color online) Spectral functions of the Hubbard model 
with the 2D square lattice DOS in the magnetic field; half-filling 
case. The inset shows the field-dependence of the magnetization. 



goes to one with increased doping. It should be noted that at 
half-filling, ReS(ti;) has a diverging slope, thus Z ^ 0. The 
system is thus not a regular Fermi liquid, but rather a singular 
Fermi liquid'*' . 
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Figure 3: (Color online) Spectral functions for the Hubbard model 
with the 2D square lattice DOS for different doping levels. The 
curves are offset vertically for clarity. The inset shows the quasi- 
particle residue as a function of the doping. 

The system exhibits unusual properties in the magnetic 
field. At half-filling, the behavior of the quasiparticle peak 



At finite filling, the behavior in the magnetic field becomes 
more in line with that of systems with non-diverging DOS, see 
Fig.|5] With increasing field, the quasiparticle peak sphts. The 
majority-spin resonance is reduced in amplitude and eventu- 
ally disappears as the lower Hubbard band becomes increas- 
ingly non-interacting-like. At the same time, the minority- 
spin resonance becomes wider and transforms into a feature 
with a sharp edge near the Fermi level. The spin-minority up- 
per Hubbard band is significantly renormalized even for rela- 
tively strong magnetic fields. Such behavior is analogous to 
that found for the Bethe lattice^. 



V. DOS WITH POWER-LAW SINGULARITY 

Motivated by the unusual features induced by the presence 
of a logarithmic (i.e. rather mild) singularity at the Fermi level 
in the DOS of the 2D square lattice, we now study the case of 
stronger power-law singularities on the example of an inverse 
square root DOS: 
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-1/2 



(17) 



The DMFT results indicate that in the presence of the interac- 
tion, the power-law singularity remains, however its exponent 
changes: the spectral functions feature a w"" singularity with 
the exponent around a sa 0.40. The energy below which this 



5 



S 2 

Q. 

1 





1.5 

1 1 

Q. 



0.5 



- 






— 0.02 




— 0.04 




— 0.08 




— 0.12 




— 0.16 




— 0.18 




\ 0.19 


-^i 1 1 1 1 1 1 






1 UAV=1.2,(n)=0.9 







a) 



-0.5 



0.5 
coAV 



1.5 



Figure 5: (Color online) Spectral functions of the Hubbard model 
with the 2D square lattice DOS in the magnetic field; finite-doping 
case. 



power-law scaling of the spectral spectral function holds de- 
pends exponentially on the value of U for small U, while for 
U > O.IW, it starts to hold essentially on the scale of bare 
parameters (i.e. U itself). In Fig. |6]3,c we plot the modified 
spectral functions^ 



(18) 



which reveal the structure of the power-law singularity. The 
decrease of F{Q) with increasing U corresponds to the pro- 
gressive naiTowing of the quasiparticle peak as the MIT is ap- 
proached. The transition is found to occur for Uc/W = 1.32 
or Uc/Wcff — 1.48, again in good agreement with the stan- 
dard result. As for other DOS functions, the Mott MIT pro- 
ceeds by the same route-"^-"*"*: with increasing U a region of re- 
duced spectral density appears, while the weight of the quasi- 
particle peak decreases until the peak disappears at the critical 

Uc. 

Interestingly, at half filling the effective quantum impurity 
model is a pseudo-gap Anderson mode l"*^^'"^^-"^^-"^^'"^^-"^^-^"'^ ' '^^ : 
the hybridization function r(a;), shown in Fig. |7h, has a 
power-law pseudogap \uj\'' with the exponent r ~ 0.38. This 
is a direct consequence of diverging spectral function and the 
imaginary part of the self-energy ImS(cj) going to at w = 0. 
Since lmS{uj) has a cusp-like cj'^ singularity with the ex- 
ponent A « 0.80, see Fig. |7j3, the self-energy goes to zero 
faster than the hybridization function (r < A), thus the system 
can be classified as a generalized Fermi liquid in the sense 
of Refs. 49 and 43. Within numerical accuracy, the spectral 
function exponent a is equal to the hybridization exponent r. 
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Figure 6: (Color online) a) Spectral functions for the Hubbard model 
with the power-law DOS in the paramagnetic phase at half- 

filling. b,c) Reseated spectral functions on linear and logarithmic 
energy scale. The arrows show the direction of increasing U. 



as expected^*^'^'. 

The low-energy expansion of the self-energy can be ex- 
pressed as^ 

E(cj) cx -\uj\^ [i + tan {X{tt/2)) sgncj] (19) 

Since the inverse spectral function goes to zero faster than 
Lu and the self-energy the self-consistency condition 

Eq. ( [Tol l reduces to 



(20) 



Furthermore, taking into account the non-interacting DOS, the 
local Green's function can be expressed as 



Gioc(t^) ^ Gq[uj - T.{uj)] 



(21) 



which reduces to G{uj) = Go [— S(cj)], since T^{uj) goes to 
zero slower than uj. If Gg [z] is the free-electron propagator in 
a system with power-law DOS It^jl^^, the low-ui expansion of 
the hybridization function is found to be 



r((jj) (X \u! 



XR 



(22) 



In our case with R ~ 1/2, we have A = 0.8 and the nu- 
merically obtained exponent r ~ 0.38 agrees well with the 
expected value of r = XR = 0.4. This implies that r < A 
for all i? e [0 : 1] for which the singularity is integrable, thus 
a strong-coupling fixed point (rather than local-moment fixed 
point) is expected in general and the system will always be a 
generalized Fermi liquid. Furthermore, the perturbation the- 
ory in U for the pseudo-gap Anderson impurity model gives a 
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relation A = 2 — 3r— , thus we obtain a result for the exponent 
of the self-energy 



A = 



l + 3i?' 



(23) 



which is expected to hold in general. For R = 1/2, this yields 
A = 4/5, which is corroborated by our numerical results. 




Figure 7: (Color online) The hybridization function and the imagi- 
nary part of the self-energy for the Hubbard model with the power- 



law DOS 



-1/2 



in the paramagnetic phase at half-filling, T = 0. 



The insets are close-ups on the low-energy region. 



The behavior in the magnetic field is very different from 
that found in previously known cases, see Fig. [8] With in- 
creasing field, the diverging quasiparticle peak transforms into 
a finite peak slightly away from the Fermi level and its width 
grows with field. At the same time, the majority-spin lower 
Hubbard band becomes increasingly non-interacting-like with 
emerging van Hove singularities in its center and sharp band 
edges. The MIT is not driven by quasi-particle mass enhance- 
ment and vanishing width of the quasiparticle band. Instead, 
it occurs when the edge of the quasiparticle band crosses the 
Fermi level. 

Upon doping, the quasiparticle peak transforms smoothly 
into the van Hove singularity in the center of the empty 
band, see Fig. |9] At finite doping, the effective hybridization 
function r{uj) no longer attains zero value at uj — Q, how- 
ever it still exhibits a sharp pseudo-gap feature with a mini- 
mum which is shifted away from the Fermi level and where 
r{uj) > 0. For small doping the minimum still appears cusp- 
like, however it becomes increasingly parabolic-like for larger 
doping. The evolution of the quasiparticle residue as a func- 
tion of the doping is shown in the inset to Fig. |9] At half- 
filling, the system is a generalized Fermi liquid, thus Z — 0. 
At finite doping, the system is a genuine Fermi liquid with a 
self-energy which behaves as ImE ^ uj"^ near the Fermi level. 
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Figure 8: (Color online) The spectral function of the Hubbard model 
with the power-law DOS |(.<j|~^'''^ in the magnetic field. The curves 
are offset vertically for clarity. The inset shows the field-dependence 
of the magnetization. 
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Figure 9: (Color online) Spectral functions for the Hubbard model 



with the power-law DOS 



-1/2 



for different doping levels. The 



curves are offset vertically for clarity. The inset shows the quasipar- 
ticle residue as a function of the doping. 



VI. CONCLUSION 

We have shown that the systems with a non-interacting den- 
sity of states which diverges at the Fermi level behave as sin- 
gular Fermi liquids, since the singularities are not washed out 
by the interactions. Their distinguishing characteristic is the 
different behavior in the magnetic field, in particular the ab- 
sence of field-induced quasiparticle localization. This has im- 
plications for the possibility of fully polarizing such systems 
with external magnetic fields of the order of the width of the 
quasiparticle band. 
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